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ABSTRACT 

We present preliminary source counts at 6.7/im and 15/im from the Pre- 
liminary Analysis of the European Large Area ISO survey, with limiting 
' flux densities of ~ 2mJy at 15/im and ~ ImJy at 6.7/xm. We separate 

^ . the stellar contribution from the extragalactic using identifications with 

APM sources made with the likelihood ratio technique. We quantify the 
completeness and reliability of our source extraction using (a) repeated 
observations over small areas, (b) cross-IDs with stars of known spec- 
tral type, (c) detections of the point spread function wings around bright 
sources, (d) comparison with independent algorithms. Flux calibration at 
15/im was performed using stellar identifications; the calibration does not 
agree with the pre-flight estimates, probably due to effects of detector 
hysteresis and photometric aperture correction. The 6.7/im extragalactic 
counts are broadly reproduced in the Pearson & Rowan-Robinson model, 
but the Franceschini et al. (1997) model underpredicts the observed source 
density by ~ 0.5 — 1 dex, though the photometry at 6.7/im is still prelim- 
inary At 15/im the extragalactic counts are in excellent agreement with 
the predictions of the Pearson & Rowan-Robinson (1996), Franceschini et 
al. (1994), Guiderdoni et al. (1997) and the evolving models of Xu et al. 
(1998), over 7 orders of magnitude in 15/zm flux density. The counts agree 
with other estimates from the ISOCAM instrument at overlapping flux 
densities (Elbaz et al. 1999), provided a consistent flux calibration is used. 
Luminosity evolution at a rate of (1 + z) 3 , incorporating mid-IR spectral 
features, provides a better fit to the 15/im differential counts than (1 + z) 4 
density evolution. No-evolution models are excluded, and implying that 
below around 10 mJy at 15/zm the source counts become dominated by 
an evolving cosmological population of dust-shrouded starbursts and/or 
© 1999 RAS act i ve galaxies. 
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1 INTRODUCTION 

The IRAS mission enjoyed huge successes, includ- 
ing the sensational discoveries of ultra- and hyper- 
luminous galaxies and of an enormous population of 
evolving starbursts. However, the survey had several 
drawbacks. For example, the bright limiting flux den- 
sities restricted the samples to low redshifts (z < 0.3) 
for all but a few ultraluminous objects. Also, only 
~ 1000 galaxies were detected at 12^tm over the whole 
sky. These deficiencies restricted the study of IR- 
luminous galaxies at all redshifts. 

The Infrared Space Observatory (ISO) offered 
~ x 1000 improvements in sensitivity in the mid-IR 
over IRAS, and the large allocations of guaranteed 
and discretionary time for deep surveys on ISO will 
greatly improve on the IRAS surveys in the mid-IR. 
For instance, ISO observations of the northern Hubble 
Deep Field (Serjeant et al. 1997, 1999, Goldschimdt et 
al. 1997, Oliver et al. 1997, Aussel et al. 1999, Desert 
et al. 1999) reached the 15/im confusion limit (~ 0.1 
mjy) over 17 square arcminutes, while the CAM-Deep 
and CAM-Shallow surveys (Elbaz et al. 1998a,b) were 
slightly less sensitive but had wider areal coverage 
(0.5 mjy over 0.3 square degrees and 0.8 mjy over 
0.41 square degrees). These have also been comple- 
mented by deep ISO photometry of selected high-z 
galaxies (e.g. Flores et al. 1999). 

The European Large Area ISO Survey (ELAIS, 
Oliver et al. 1999 (paper I), Rowan- Robinson et al. 
1998) was the largest open time project on ISO, com- 
plementing the deep ISO samples by surveying ~ 12 
square degrees to a depth of ~ 2 mjy at 15/im and 
< 100 mjy at 90/im. Around half the area was also 
mapped at 6.7/zm to ~ 1 mjy. Three fields in the 
Northern hemisphere (Nl, N2, N3) collectively com- 
prised around two-thirds of the 15/im areal coverage, 
with the remaining area taken by the Southern SI 
field and several small areas in both hemispheres. The 
ambitious cosmological aims include tracing the ex- 
tinguished star formation history of the Universe to 
z ~ 1 — 2, orientation-independent selection of dust- 
shrouded quasars, and the potential discovery of hy- 
perluminous galaxies (with comparable intrinsic lu- 
minosities to IRAS FSC 10214+4724) out to redshifts 
z S 5. A more detailed discussion of the diverse scien- 
tific aims of ELAIS, the selection of areas and observa- 
tional parameters can be found in the ELAIS survey 
paper (Oliver et al. 1999); in summary, the survey 
areas were selected to have low galactic cirrus emis- 
sion, high visibility by ISO, high ecliptic latitude and 
avoiding 12/im IRAS sources brighter than 0.6 Jy. In 
another companion paper, Efstathiou et al. 1999, we 
discuss the 90/j.m source counts from the Preliminary 
Analysis of the ELAIS ISOPHOT data, and in Crock- 
ett et al. (1999) we discuss the stellar mid-infrared 
source counts. The ELAIS areas have also been the 
subject of intensive multi-wavelength follow up, sum- 
marised to date in Oliver et al. (1999) and presented 
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in detail in other papers (e.g. Ciliegi et al. 1999, Grup- 
pioni et al. 1999). Here we present the completeness, 
reliability and extragalactic source counts from our 
initial Preliminary 6.7/xm and 15/j.m ISOCAM cata- 
logues. A future paper will present the Final Analy- 
sis products from the ISOCAM ELAIS data, which 
is expected to improve on the Preliminary Analysis 
presented here. 

This paper is structured as follows. In section |^ 
we describe the Preliminary Analysis CAM pip eline, 
explaining the artefacts in the data (section 2.1), and 
the pipeline algorithm (section 2.2). The results from 
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the Preliminary Analysis catalogue are presented in 
section ^[ Our various completeness and reliability es- 
timates are discussed in section and the segrega- 
tion of extragalactic from stellar sources in section 3. J . 
Section 1.4 presents the source counts in both wave- 
bands. These results are compared with source count 
models and previous results in section ^, where we 
also discuss the implications for the evolution of star 
forming galaxies and on the star formation history of 
the Universe. 



2 ELAIS CAM PRELIMINARY 
ANALYSIS 

2.1 Data quality 

The ELAIS CAM survey proper was conducted in 
raster mode (astronomical template CAM01), with 
the LW-2 (6.7/im) and LW-3 (15/xm) filters. Details 
of the CAM01 Astronomical Observation Template 
(AOT) can be found in paper I. The CAM detector is 
stepped across the sky in a grid pattern, with roughly 
half-detector-width steps in one direction and roughly 
whole detector widths in the other, covering approx- 
imately half a square degree per raster. This pattern 
leads to a redundancy of at least x 2 over most of the 
area surveyed. At each raster pointing (i.e. each grid 
position of the raster) the 32 x 32 CAM detector is 
read out several times. 

Like the ISO-HDF North data (Serjeant et al. 
1997, 1999, Aussel et al. 1999, Desert et al. 1999) 
the ELAIS CAM data contains many problematic 
artefacts. Because of the frequent and complicated 
glitches, we do not take the approach of reconstruct- 
ing a sky-map and searching for sources in these maps. 
Rather, we look for the characteristic signatures of 
sources and glitches in the time histories of individ- 
ual pixels. See Starck et al. (1998) or Aussel et al. 
(1999) for more details. 

The CAM detector also exhibits hysteresis. 
Source fluxes are initially around a factor of 2 fainter 
in instrumental units than the stabilised (i.e. asymp- 
totic) value. Our survey strategy ensures that sources 
almost always have corroborating sightings in sep- 
arate pixels. This permits a filter to remove glitch 
events from candidate sources. 
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Figure 2. Effective areal coverage in the LW-3 15/nn 
band. The maximum coverage, which excludes the first 
and last pointings of each raster and regions without re- 
dundancy (removing 7% from each raster), is O ma x = 10.0 
square degrees. The curves asymptote to < 1 because the 
bright end of the PSF wing test is not well constrained, 
having only a few weak limits. Also the asymptotic value 
at low fluxes is not included (and is also ill-constrained) 
since this is presumably glitch confirmation. 



2.2 Preliminary source extraction pipeline 

2.2.1 Preprocessing and deglitching 

The available CAM data reduction software under- 
went several substantial improvements over the life- 
time of ISO, as the knowledge of the detector char- 
acteristics improved. However, from the outset we 
needed a method of preliminary data analysis, to feed 
for example immediate follow-up projects. Such a Pre- 
liminary Analysis pipeline may not of course repre- 
sent best practice at the end of the ISO mission, but 
should at least provide reasonably complete and reli- 
able preliminary source list. It was decided that the 
CAM Preliminary Analysis data reduction should be 
as uniform as possible, which required that the Pre- 
liminary Analysis pipeline be fixed at an early stage. 
Accordingly our adopted pipeline could not incorpo- 
rate the accurate field distortion corrections (Abergel 
et al. 1998), which were not established at the start of 
the mission, nor the analytic models for the cosmic ray 
transients and detector hysteresis which were devel- 
oped in the course of the ISO lifetime (e.g. Lari 1999, 
Abergel et al. 1998). Nevertheless, such improvements 
will be incorporated in future ELAIS Final Analysis 
products. 

The data reduction was performed using the In- 
teractive Data Language (IDL) software. The ini- 
tial steps of the CAM Preliminary Analysis pipeline 
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Figure 1. Panels show the 15/im (left panels) and 6.7/xm (right panels) flux calibration data. Top panels show the 
distribution of the fraction of flux in the brightest (6") pixel, out of two randomly placed observations of a point source. 
To calculate this, we use an evenly-spaced grid of theoretical PSF models from Aussel (priv. comm.). The distribution is 
skewed to higher fluxes compared to the expected peak flux distribution of a single observation, and has smaller variance. 
The central panels show the detected fraction of individual pixels in the PSF wings of bright stellar sources. The central 
pixels in each star have been excluded from these plots. The lines show the best fit parametric model using equation [lj 
The lower panels show comparisons of the detected fluxes in individual pixels in bright point source PSF wings, with the 
predictions based on the PSF model. Again, the central pixels have been excluded from these plots. 
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Figure 3. Effective areal coverage in the LW-2 6.7^tm 
band. The maximum coverage, which exciudes the first 
and last pointings of each raster and regions without re- 
dundancy (removing 7% from each raster), is n max = 6.51 
square degrees. 



Figure 4. Completeness at 15/^m estimated from the 
predicted fluxes of stars in the main ELAIS regions. For 
more details see Crockett et al. 1999. 



are straightforward. The edited raw data supplied by 
ESA were converted to IDL structures using the CAM 
Interactive Analysis (CIA, April 1996 version), and 
converted from ADUs (analogue-to-digital units) to 
ADUs per second. The default dark frame was sub- 
tracted from each exposure. To estimated the noise 
level in each pixel, we performed an iterative Gaussian 
fit to the histogram of readout values for each pixel. 
Cosmic ray spikes were then identified by > 4a rises 
followed (one or two readouts later) by > 4<r falls. 
A similar algorithm was used to identify occasional 
readout troughs. The readout histograms were then 
re- fit and an empirical sky flat field was obtained; the 
default ESA-supplied flat field was found to give very 
unsatisfactory results. An attempt was made to model 
the initial detector stabilisation using the IAS model 
within the CAM Interactive Analysis (CIA) software 
package. 

2.2.2 Background estimation and Source Detection 

Unlike the detectors on the IRAS satellite (e.g. Neuge- 
bauer et al. 1984), the transients in the ISOCAM 
detector make the background levels in each pixel 
vary strongly and discontinuously with time. The ap- 
proach adopted in the Preliminary Analysis was to 
estimate the background level in a given pixel and 
pointing from linear fits (in time) to the readouts in 
previous and subsequent pointings, then identify can- 
didate sources from > 3<r features above the back- 
ground in the pixel readout timelines. The source ex- 



traction pipeline is therefore spatially 1 dimensional. 
This simple approach avoids any explicit parameteri- 
sations of the transient profiles, which were not avail- 
able at the time, and allows a local error estimate of 
the background level. This error was taken to be the 
formal error on the fit, but did not incorporate the 
instrumental noise at each data point. Since cosmic 
rays were frequently observed to cause discontinuities, 
the range for the linear fits extended no further than 
the nearest (in time) cosmic ray, and in any case not 
longer than 3 pointings. If an acceptable \ 2 was n °t 
obtained in the linear fit, the range for the fit was 
decreased; if a good fit was still not obtainable, an 
average of the 10 nearest readouts was taken. Sources 
could then be identified from their excess above the 
extrapolated backgrounds. Since sources also create 
discontinuities, the data stream was iteratively re- 
fit treating > 3.5<r sources in the same way as cos- 
mic rays above. An initial list of brightest sources 
was obtained before the iterative fitting, by searching 
for their discontinuous rise at the start of a point- 
ing, and discontinuous fall at the end. Where sources 
were found (whether initially or in the iterative fit- 
ting) the background level in later iterations was ex- 
trapolated from only the previous pointings, ignoring 
the subsequent pointings, since the hysteresis after a 
source would otherwise lead to an overestimation of 
the background level. Sources were not extracted from 
the first or final pointings of the rasters, owing to lim- 
itations in the background fitting routines at the time 
the pipeline was frozen. Note that this iterative source 
extraction does not distinguish genuine sources from 
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Figure 5. Flux calibration at 15/mi estimated from the 
predicted fluxes of stars in the main ELAIS regions. The 
observed fluxes assume a one-to-one conversion between 
ADU/gain/second and mJy. Open points are from Hip- 
parcos, and solid points from Simbad. As discussed in 
Crockett et al. 1999, the most reliable spectral typing and 
hence mid-IR predictions are available for Simbad stars. 
The straight lines show a 1 : 1 and 2 : 1 ratio between ob- 
served and predicted counts. Also shown is the ±1<t limits 
on the mean observed:predicted ratio as a function of flux, 
calculated within ±0.25 dex of each flux (i.e., an 0.5 dex 
boxcar smoothing). There are perhaps hints from this fig- 
ure that the flux calibration may be a weak function of 
flux, though in this paper we assume a flux-indcpcndcnt 
scaling. For more details see Crockett et al. 1999 and Mis- 
oulis et al. 1998. 



"fader" transient events, so a further filter for source 
candidates is still required. 



2.2.3 Source Corroboration 

The source detections in the pixel readout histories 
were spatially merged in each pointing using a con- 
nected pixel algorithm. (Note that there was no mini- 
mum number of pixel detections, unlike e.g. connected 
pixel source detection on an oversampled CCD im- 
age.) We then used the ~ x2 redundancy in the 
CAM rasters to search for corroborating observations 
of each source candidate. Genuine sources should be 
present in both observations, but glitch events should 
not be confirmed except by chance. The adopted 
search radius of 2 pixels, while large enough to safely 
encompass the (then uncertain) field distortion, nev- 
ertheless led to a large number of spurious detections, 
with the majority of source candidates at 15/im due 
to glitch corroboration. Each candidate corroborated 
source was therefore examined by eye independently 



by at least two observers, who assigned quality flags 
of 1 — 4, where 1 refers to a "definite source," 2 is a 
"probable source," 3 "probably not" and 4 "definitely 
not." Each raster in the LW-2 (6.7/im) filter yielded 
typically ~ 100 events in total of which around one 
half were classified as probable sources by at least one 
observer. In LW-3 (15/im) both the number and frac- 
tion of spurious events was much higher: there were 
typically ~ 30 strong source candidates, and a further 
30—50 sources where the classification was ambiguous 
or debatable, with typically around 200 — 500 spurious 
events. Note that the Preliminary Analysis algorithm 
will necessarily miss sources with (a) only one obser- 
vation or (b) corroborating observations in only the 
first or last pointings, so the effective area is slightly 
less than the nominal ~ 12 square degrees. 



2.2.4 Astrometric corrections 

After the Preliminary Analysis reduction was com- 
plete, we improved the astrometry by incorporating 
the latest field distortion correction into the 15/im 
source catalogues. Several sources with strong tran- 
sient events nearby had their centroids strongly af- 
fected by the glitches. We therefore adopted a sim- 
ple strategy for our Preliminary Analysis astrome- 
try and flux calibration: the flux and (distortion- 
corrected) position of a source are taken from those 
of the brightest single-pixel detection of that source, 
excluding transients. We found this to be superior to 
(eg) masking nearby transients by eye then recalcu- 
lating the centroids of the eyeball-accepted sources, 
particularly if the PSF wings lie on the detector but 
the source itself is just outside. Our adopted algo- 
rithm should yield astrometry accurate to ±3" in the 
absence of any other systematic errors. Two such sys- 
tematics were expected in our data: firstly, errors in 
the position of the lens introduce a random astro- 
metric offset to each raster of order one pixel; sec- 
ondly, any errors in the calculation of the pointing po- 
sition by CIA would offset any sources in that point- 
ing. By examining the offsets with the likelihood-ratio 
identifications we can determine the lens offset em- 



pirically (section 3.3). However, several rasters were 
found to have bimodal distributions of ISO-optical 
offsets, due to some unknown error in the CIA-derived 
astrometry in at least part of the raster. We there- 
fore rederived the pointing astrometry using the ESA- 
supplied IIPH.FIT astrometry file, using the median 
coordinate positions in the duration of the pointing. 
This was found to remove the bimodality. 



2.2.5 Aperture Corrections to Photometry 

As discussed above, our source extraction method in- 
volves looking for characteristic time signatures in in- 
dividual pixels. Without (at the time) a reliable and 
exact model for the glitch events, nor a reliable glitch 
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event identification, we found that aperture photom- 
etry around our source positions was often seriously 
affected by nearby glitches. Instead of aperture pho- 
tometry, we simply took the brightest flux of the pix- 
els detecting the source, excluding (by eye) those pix- 
els affected by glitches. Clearly, some aperture cor- 
rection is needed to correct these peak pixel fluxes to 
total fluxes. 

We can quantify these aperture corrections using 
a PSF model. In figure |l| we show the predicted flux 
in the brightest pixel of two randomly positioned ob- 
servations of a point source (recall that at least two 
observations of a source are required for it to pass the 
Preliminary Analysis selection). At 15/xm, the bright- 
est pixel has a flux of ~ 0.4 ±0.1 times the total 
flux of that source. At 6.7/im the histogram is more 
sharply peaked, since the PSF is undersampled; the 
peak flux is always less than 0.69 times the total at 
6.7/im, but is greater than 0.5 (0.6) of the total in 
> 90% (> 75%) of occasions. We therefore applied 
global aperture corrections of 2.36 at 15/im and 1.54 
at 6.7/im to our peak fluxes. 



3 RESULTS 

3.1 Eyeballing results 

Our eyeballing results imply our catalogue is highly 
reliable to at least 3 mjy at 15pm. In the ELAIS 
areas considered in this paper, there were 715 15/im 
sources accepted by two obse rver s, of which 510 had 
APM identifications (section 3.3); of the 816 singly- 



accepted sources at 15/tm, 212 had APM identifica- 
tions. The singly-accepted sources are heavily skewed 
to faint fluxes: 90% are fainter than 2.2 mjy. If we 
choose to regard all sources accepted by only one ob- 
server as false positives, and all sources accepted by 
two as true positives, then we obtain ~ 80% reliabil- 
ity for sources accepted by any observer at ~ 3 mjy 
and ~ 95% at 5 mjy. These are one of the most pes- 
simistic assumptions we could make; if alternatively 
we merge the doubly-accepted sources with the op- 
tically identified singly-accepted sources, the reliabil- 
ity (fraction doubly-accepted) of the sources accepted 
by any observer rises to ~ 95% at ~ 3 mjy and 
increases at brighter fluxes. The blank-field singly- 
accepted sources are heavily skewed to faint fluxes, so 
we can neglect their effect on the counts above 3 mjy 
and combine the doubly-accepted sources with the op- 
tically identified singly-accepted for the purposes of 
the source counts. However this is not to say that the 
fainter singly-accepted sources necessarily have lower 
reliability, since the genuine fainter sources may have 
fainter optical counterparts. The situation is similar 
at 6.7/im. Of the singly-accepted sources 95% have 
fluxes less than 1.5 mjy. 



3.2 Completeness, reliability and flux 
calibration 

3.2.1 Repeated observations 

There are several potential approaches to estimating 
the Preliminary Analysis completeness and reliability, 
the most obvious of which is repeated observations 
over small areas. Oliver et al. (1999) present 10 repeat 
observations of a small ELAIS raster, covering at least 
six known sources by design. The source extractions 
in the individual reductions confirmed our result that 
the source extraction is highly complete and reliable 
at flux densities above ~ 2.5mJy at 15/im. 



3.2.2 PSF-wing test 

One robust test of our temporal source extraction is 
to examine the point spread function wings of bright 
sources. Like the cross-correlation with bright stars, 
this relies on the knowledge that a given pixel is il- 
luminated by a known source flux, but has the ad- 
vantage that the comparison can be taken to arbi- 
trarily faint flux levels. We selected the 30 brightest 
15/im sources with stellar optical identifications, and 
selected a model theoretical Point Spread Function 
(PSF; Aussel, priv. comm.) with the smallest RMS 
difference in the central 3x3 PSF pixels. We nor- 
malised the PSF to the source flux using the mean 
observed flux in the central pixel. Using this model 
we predicted the mean flux in each pixel of the PSF 
wings, and hence determined the single-pixel detec- 
tion efficiency / of the temporal source extraction. 
The results are shown in figure We can also compare 
the predicted PSF fluxes with the fluxes extracted b 
the Preliminary Analysis algorithm, shown in figure 
The relation is encouragingly linear, though the scat- 
ter is larger than expected (by ~ x2 at 15/im) based 
on the quoted errors in the Preliminary Analysis sky 
background fits. This is perhaps not surprising, since 
the background fitting algorithm does not make use 
of the detector noise in the fit, so will tend to un- 
derestimate the background level uncertainty; also, 
non-white noise features may prevent the noise scal- 
ing as y/N. Oddly, the discrepancy in the scatter is 
largest at brightest fluxes. Plausible explanations in- 
clude slight errors in the source centroids or in the 
theoretical PSF shape, both of which sensitively af- 
fect the brightest flux predictions. Analogous results 
for the 30 brightest 6.7/im stars are also shown in fig- 
ure []]. The undersampling of the PSF at 6.7/im makes 
the scatter harder to interpret: more of the flux is con- 
tained in the central pixel making it less sensitive to 
errors in the assumed PSF shape, but is more affected 
by the much larger uncertainty in the centroid. 

There are several caveats which apply to the high 
apparent completeness in figure [l]. The chance detec- 
tions of "faders" will tend to increase this estimate, 
but since both the 6.7/im and 15/im completeness ap- 
pear to asymptote to ~ 0.015 (i.e. probability of de- 
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tection of nearby spurious events is 1.5%) this appears 
not to be a serious problem. This is also only appli- 
cable to single-pixel detections, whereas in fact the 
CAM detector Nyquist samples the PSF at 15/^m. 
Indeed at 15pm the fraction of the flux in the bright- 
est pixel rarely exceeds 0.5 though the same does not 
apply at 6.7pm. The completeness also is not the de- 
tected source fraction F but rather F 2 , since we re- 
quire at least 2 independent detections for a source to 
be accepted by our algorithm. There is also a slight 
bias in that the brightest detected sources are not 
typically observed when the detector is suffering the 
strongest transients, because the sources would not 
otherwise be detected. Finally, any incompleteness 
caused by the eyeballing, or by any areal coverage lost 
to e.g. cosmic ray impacts, would not be included in 
these estimates as they stand. 

We can therefore estimate the Preliminary Anal- 
ysis completeness from the PSF-wing test in the fol- 
lowing way. We can convolve the single-pixel detec- 
tion rate in figure |l| with the peak flux distribution 
expected from the PSF models to obtain the source 
sensitivity F for single observations. The Preliminary 
Analysis completeness, before eyeballing, will then be 
proportional to F 2 , assuming the PSF wings them- 
selves are representative of the data as a whole. 

We fit the detected fractions in figure [I] with tanh 
functions, i.e. 



f{S] 



{fn 



/min) 



tanh(alog 10 (S/S ) + l)+/min(l) 



where / m i n and / max represent the asymptotic lim- 
its at faint and bright fluxes respectively. We de- 
fine the single-pixel source detection fraction to be 
f'(S) = f(S) - /nun, and use a grid of PSF models 
(each normalised to 1) spanning the possible range 
of centroid positions to estimate the single-pointing 
source detection: 



x S, 



pcak,i ) 



(2) 



where Speak, i is the peak flux in the ith PSF. This 
assumes that if a source is not detected in the peak 
pixel, it will not be detected in any pixel. The max- 
imum areal coverage of the Preliminary Analysis, ie 
excluding the first and last pointings of each raster, 
and regions with no redundancy, is f2 max = 10.0 
square degrees at 15pm, and fi max = 6.51 square de- 
grees at 6.7 pm. From this we obtain the Preliminary 
Analysis areal coverage as a function of flux: 



Slprcl iminary Analysis (5) = n m ^F 2 (S) 



(3) 



The final areal coverage from the PSF-wing test is 
plotted in figures ^ and |l| 

3.2.3 Comparison with IAS and CEA pipelines 

As a final check of the completeness of our Prelim- 
inary Analysis catalogue, we compared our source 



extraction in the repeat obser vation regions (c.f. 
Oliver et al. 1999 and section 3.2.1) with extrac- 



tions made with the CAM-Deep pipeline developed 
at the Commisariat a L'Energie Atomique, Saclay 
(CEA; Baker, priv. comm.) and a pipeline based on 
the "Triple Beamswitch" method developed at Insti- 
tut d'Astrophysique Spatiale (IAS; Clements, priv. 
comm; Desert et al. 1998) . Of our six robust sources 



(section 3.2.1), CEA and IAS both identify three (the 
same three), with fluxes brighter than ~ 3mJy. This 
appears to be mainly because the Preliminary Analy- 
sis pipeline extracts lower signal-to-noise sources, but 
supplements with greater manual eyeballing. Never- 
theless, this confirms that we are not significantly un- 
derestimating the surface density of sources brighter 
than 3mJy. 



3.2.4 Flux calibration and Cross- correlation with 
bright stars 

Another useful test of the completeness is to search 
for detections at the locations of bright stars (see sec- 
tion 



3.3 for details of the optical identification algo- 



rithm). If the spectral types of the stars are known, 
one can predict their mid-infrared fluxes. An excep- 
tion would be dust-shell stars, but these are expected 
to be rare in the survey. In Crockett et al (1999) and 
Crockett (1999) the sources are cross-correlated with 
stars from the Simbad and Hipparcos databases. All 
the 22 stars with predicted fluxes brighter than 3mJy 
at 15pm were detected in the Preliminary Analysis, 
and all but two of the 20 stars with predicted 6.7pm 
fluxes above lmjy appeared in the Preliminary Anal- 
ysis. To assess the level of random associations, we 
randomised the stellar positions within the ELAIS 
fields and repeated the cross-association, and found 
none. The 15pm completeness is shown in figure ^. 
Note that this is an extremely robust estimator of 
the completeness, since the stellar flux prediction al- 
gorithm reproduces IRAS fluxes well, and the stars 
cannot be argued to lie on atypical regions of the de- 
tector. 

However, the flux calibration implied by these as- 
sociations is around a factor of 2 discrepant with the 
expectation that ADU/gain/second ~ mjy at 15pm 
(see figure |^). Across the entire range in flux, the 
±1<7 limits on the log of the calibration ratio are 
0.246±0.050, ie the ratio is 1.75±^;§|. As shown in fig- 
ure H there are hints that this calibration is a function 
of flux, with a calibration ratio of 2 preferred at faint 
fluxes (where indeed most of our sources lie). (The 
ISOCAM Observer's Manual recommends a conver- 
sion of approximately 2 ADU/gain/second to 1 mjy 
at 15pm for fully stabilised sources, and after correct- 
ing for the loss of flux due to lack of stabilisation (e.g. 

To 



2.1) becomes around a 1 : 1 conversion. 



section 

predict the fluxes we followed the procedure of Crock- 
ett et al. (1999) and Misoulis et al. (1998), incorpo- 
rating the passband profiles (for more details we refer 
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the reader to these papers; the algorithm accurately 
reproduces stellar mid-IR IRAS fluxes so there are un- 
likely to be significant systematics in the CAM flux 
predictions). It is not immediately clear what might 
cause such a discrepancy, though there are several 
possibilities, such as the uncertainties in the PSF for 
the aperture c orre ction, and the assumed ~ x2 loss 
in flux (section 2.1) from the lack of an upward source 



LR(f,x,y) oc 



n{I) 



(4) 



stabilisation correction. (Note that differences in the 
PSF due to the differing spectral slopes of the stars 
are too small to account for the discrepancy: e.g. Aus- 
sel et al. (1999) find it only to be a ~ 10% effect.) For 
the purposes of the source counts we will adopt the 
mjy : ADU/g/s = 1:2 stellar calibration implied in 
figure where the ADUs are not corrected for losses 
due to lack of stabilisation. 

At 6.7pm a lower bound of 95% can be made on 
the completeness at fluxes > 10 mjy, but the uncer- 
tain aperture corrections make applying the stellar 
flux calibration more difficult. As a preliminary mea- 
sure we therefore simply take the pre-flight ISOCAM 
Observer's Manual cal ibra tion at 6.7pm, corrected by 
a factor ~ 2 (section 2.1) to account for the loss in 
flux from lack of stabilisation. 

In summary, our various completeness estimates 
yield a > 1 mjy limit to the completeness at 15pm, 
and > 0.5mJy at 6.7pm (figures [| and |^) using the 
ISOCAM observer's manual flux calibration corrected 
by a fact or of 2 to account for stabilisation loss (sec- 
tion 2.1). However our stellar cross-correlation sug- 
gests we have underestimated the fluxes by a factor 
of ~ 2 at 15pm (figure ^) so the 15pm completeness 
quoted should by revised upwards to > 2 mjy. In 
all subsequent discussion, the 15pm ELAIS fluxes are 
assumed to obey this stellar flux calibration. 



3.3 Optical identifications 

In this section we summarise the optical identifica- 
tion algorithm used for the Preliminary Analysis cat- 
alogue, for APM stars and galaxies (McMahon & Ir- 
win 1992). We adopted the likelihood ratio proce- 
dure of Sutherland & Saunders (1992) to associate 
our Preliminary Analysis sources with known objects. 
The surface density of catalogue objects as a function 
of magnitude is incorporated into the likelihood of 
each identification of the Preliminary Analysis sources 
with the catalogue objects. Following Mann et al. 
(1997), we define the likelihood ratio to be the ratio 
of the probability of detecting a genuine counterpart 
to the source with the position and flux of the cata- 
logue object, to the probability of such an association 
occurring by chance given the positional errors. For a 
catalogue surface density n(f) (where / is the flux), a 
positional uncertainty e(x, y) (where x and y are e.g. 
Cartesian coordinates) and an a priori flux distribu- 
tion of IDs given by q(f), the likelihood ratio is given 
by 



Using this expression and, assuming q(f) to be con- 
stant as a function of optical magnitude, we found, 
for each Preliminary Analysis source, the object in the 
APM catalogue giving the highest likelihood ratio. We 
normalised the likelihood ratios by finding the max- 
imum likelihood-ratio associations around random 
source positions, to yield probabilities for each Pre- 
liminary Analysis identification, and accepted iden- 
tifications with random probabilities of less than 0.3. 
Slight errors in the lens positioning lead to systematic 
astrometric shifts in each raster, of order a few arc- 
seconds. To correct for this, the high-likelihood iden- 
tifications in each raster were used to determine any 
systematic astrometric offset. The identifications were 
subsequently rederived. There were not enough reli- 
able optical identifications to obtain a robust estimate 
of the lens offset in the smaller ELAIS areas, so the 
analysis was restricted to the main ELAIS areas of 
Nl, N2, N3 and SI. Further discussion of the identi- 
fications is deferred to later papers in this series. 



3.4 Source counts 

Usi ng the completeness and reliability from section 



3.2, we can extract the extragalactic source counts 
at 6.7pm and 15pm. For the purposes of the counts 
we included all sources accepted by (at least) two 
observers. From this list, we exclude 15/im stellar 
identifications brighter than approximately B — 16.5 
(assuming a monotonic stellar plate saturation cor- 
rection converting B = 18 galaxy magnitudes to 
B < 16.5 stars), but include fainter stellar IDs since 
these are expected to be predominantly AGN. All 
stellar IDs at 6.7pm were eliminated from the extra- 
galactic source counts. It is highly unlikely for stars 
with B > 16.5 to be detected at 15pm to these lim- 
its (B — 16.5 is equivalent to Sb = 1 mjy; see also 
Crockett et al. 1999 and Crockett 1999), with the ex- 
ception of rare dust-shell stars. Note also that all the 
15pm (6.7pm) sources accepted by only one observer 
are fainter than 2.2mJy (2.5mJy), so we can be highly 
confident of the reliability of sources brighter than 
these limits. 

The segregation of extragalactic from stellar 
counts is robust at 15pm, but the large fraction of 
stellar IDs at 6.7pm make it possible that some faint 
stars have been included in the extragalactic counts at 
this wavelength: of the 794 doubly-accepted sources at 
6.7pm, only 79 have stellar APM classifications faint 
enough to be accepted in the extragalactic counts. We 
eyeballed the DSS images of every optical identifica- 
tion of the 6.7pm and 15pm sources. We found the by- 
eye classifications to agree with the APM in nearly all 
cases. The resultant extragalactic counts are virtually 
indistinguishable from the automated segregations at 
both 6.7pm and 15pm. 
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Figure 6. Integral ELAIS extragalactic source counts at 15^tm. Flux densities are quoted in Jy. Shaded regions show the 
ranges spanned by ±1<t uncertainties. Also shown are the source counts and P(D) analysis from the Hubble Deep Field 
(Oliver et al. 1997, Aussel et al. 1999). The IRAS counts are estimated from the 12/xm counts as discussed in the text. Top 
Left panel shows the Pearson & Rowan-Robinson model. Dotted line shows normal galaxies, dashed line shows Seyfert I 
galaxies, dot-dashed are Seyfert lis, and starbursts are dash-dot-dot-dot. The total is shown as a full line. Top Right panel 
shows the Franceschini et al. (1994) source count model overplotted. The model spiral contribution is shown as a dotted 
line, ellipticals as a dashed line, SO as a dash-dot line, starbursts as a dash-dot-dot-dot line and AGN as a long dashed 
line. The total population model is shown as a full line. Also plotted are the Guiderdoni et al. (1997) models A and E, 
as small filled crosses and small open circles respectively. Bottom Left panel shows the evolving Xu et al. (1998) models, 
renormalised by X 1.8 to match the IRAS counts as in figure ^| The full line shows (1 + 2) 3 luminosity evolution with K- 
corrections derived from star burst models with mid-IR features, and the dashed line shows the same luminosity evolution 
without the mid-IR features. Density evolution of (1 + z) 4 with mid-IR features is plotted as a dotted line, and without 
mid-IR features as a dash-dot-dot-dot line. Bottom Right panel shows all available no-evolution models. Franceschini et 
al. (in prep.) is overplotted as a dash-dot-dot-dot line. The Xu et al. (1998) models are shown with and without the MIR 
spectral features (dash-dot and dotted respectively). All three no-evolution models have been renormalised to match the 
IRAS counts, by a factor of 0.8 for the Franceschini et al. models, and 1.8 in the case of the Xu et al. (1998) models. 
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Figure 7. Preliminary integral ELAIS extragalactic source counts at 6.7fim, using the flux calibration discussed in the 
text. Flux densities are quoted in Jy. Left panel shows the Pearson & Rowan-Robinson model overplotted, and right panel 
shows the Franceschini et al. (1997) model. Symbols as in the respective panels of figure pi 



The extragalactic source counts are plotted in 
figure ^, using the by-eye star-galaxy separation and 
the stellar flux calibration. The 15^im Lockman Hole 
ISOCAM survey data will be discussed in Elbaz et 
al. 1999. The counts from this survey at around the 
~ 5mJy level, which overlap with the ELAfS counts, 
appear significantly lower than those of ELAIS. How- 
ever this is entirely attributable to the ~ x2 differ- 
ences in assumed flux calibration. The counts are in 
excellent agreement in instrumental units or in mjy 
with the same flux calibration assumption; alterna- 
tively, a ~ 30% reduction of the ELAIS Preliminary 
Analysis 15/im flux calibration factor would also bring 
the counts into formal la agreement while remaining 
consistent with figure [lj. However this may require a 
commensurate reduction in the IRAS counts. Such a 
reduction has been argued to be necessary by Elbaz et 
al. (1999) in order to account for large-scale-structure 
effects in the Rush et al. sample. 

Also plotted are the 12/xm source counts cal- 
culated by Oliver et al. (1997) from the Rush et 
al. (1993) sample, using the QMW IRAS Galaxy 
Mask (Rowan- Robinson et al. 1991), and transposed 
to 15/xm assuming a population mix matching the 
Pearson & Rowan-Robinson (1996) predictions. The 
faint source counts from the ISO-HDF North survey 
(Oliver et al. 1997) are also shown in the counts fig- 
ures. The ELAIS extragalactic Preliminary Analysis 
counts at 15/im are consistent with an interpolation 
between the ISO-HDF North and Rush et al. (1993) 
data sets, and at 6.7/im with reasonable extrapola- 



tions from the ISO-HDF North. The source density 
at 15//m is also in good agreement with that obtained 
at 12/xm by Clements et al. (1999), though the differ- 
ing K-corrections make it not immediately clear that 
the counts must necessarily agree (e.g. Xu et al. 1998, 
Serjeant 1999). Note that the 6.7/xm counts have sig- 
nificant photometric errors independent of flux, due 
to the undersampling of the PSF (figure [l]). For a 
constant source count slope the shape of the counts is 
unaffected by flux-independent errors, so we can re- 
gard the 6.7/im counts as subject to a potential sys- 
tematic error in the form of a horizontal shift. Such 
a systematic is less than or of order a factor 2 in 
flux. We overplot the model predictions from Pear- 
son & Rowan- Robinson (1996) in figures [| and 0, as 
well as the model predictions of Franceschini et al. 
(1994, 1997). Also overplotted are the Guiderdoni et 
al. (1997) 15/im model counts and the evolving mod- 
els from Xu et al. (1998), with and without mid-IR 
spectral features. Note that all the Xu et al. models 
have been renormalised to match the IRAS counts. 
Figure |^ also compares the observations to a variety 
of non-evolving models. Apart from the renormalisa- 
tion to the IRAS Rush et al. (1993) x QIGC counts, 
these are the same non-evolving models as discussed 
in Elbaz et al. 1998. 
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Figure 8. Left panel shows Euclidean-normalised differential source counts at 15^tm. Flux densities are quoted in Jy. 
Triangles show the Rush et al. (1993) sample, and rilled circles show ELAIS. The open circle is taken from the ISO-HDF 
sample (Oliver et al. 1997). The filled line is the no-evolution model of Xu et al. (1998) with mid-IR spectral features 
incorporated. The dotted and dashed lines respectively show the Xu et al. (1998) predictions for (1 + z) 3 pure luminosity 
evolution and (1 + z) 4 pure density evolution, both including mid-IR spectral features. All three Xu et al. (1998) models 
have been renormalised by xl.8 to match the IRAS counts. Horizontal error bars (where shown) indicate the width of 
the flux bins for the fainter samples, and vertical error bars the ±lcr Poisson errors on the counts. The right panel shows 
Euclidean-normalised differential source counts at 6.7/im. Flux densities are quoted in Jy. Filled circles show ELAIS, and 
the open circle shows the ISO-HDF sources (Oliver et al. 1997) with optical IDs. 



4 DISCUSSION 

The experimental agreement with the evolving mod- 
els of Pearson & Rowan-Robinson (1996), Frances- 
chini et al. (1994), Guiderdoni et al. (1997) and Xu 
et al. (1998) at 15/xm over seven orders of magnitude 
in flux density is striking. The starbursts in the Pear- 
son & Rowan-Robinson (1996) models are normalised 
to the 60/xm IRAS counts, but not explicitly to the 
12^im counts. The slight overprediction of the IRAS 
counts is also present in the ELAIS counts, until a 
slight upturn at around lOmJy (which is not an effect 
of incompleteness or low reliability) departing from 
the Euclidean slope brings the data into closer agree- 
ment with the model. Of the four Xu et al. (1998) 
models, the (1 + z) 3 luminosity evolution models have 
the stronger upturn, reproducing the observed counts 
slightly more well than the alternative (1 + z) 4 den- 
sity evolution. This is clearer still in the 15/mi differ- 
ential counts (figure ^) , where luminosity evolution is 
shown to make a much better fit to the ELAIS counts, 
though an even larger excess is suggested by the ISO- 
HDF counts. Note that we renormalised the Xu et al. 
(1998) predictions by x 1.8 to match the IRAS counts. 

The source count models present rather differ- 
ent predictions at 6.7^im. The Pearson & Rowan- 



Robinson (1996) model has only a slight overestimate, 
accountable for instance to the flux calibration uncer- 
tainties. However, the underprediction in the Frances- 
chini et al. (1997) model is over an order of magni- 
tude, probably too large to be an artefact of our albeit 
uncertain flux calibration at this wavelength. This dis- 
crepancy is most likely to be due to deficiencies in the 
assumed spectral energy distributions, which are not 
well-constrained in this wavelength range (e.g. Ser- 
jeant 1999), rather than due to incorrect assumptions 
about the evolution or population mixes. 

The loss of the Wide Field Infrared Explorer 
(WIRE) satellite was a serious blow to infrared ex- 
tragalactic astronomy. In the hope or expectation of a 
new WIRE mission, we can compare our source counts 
with the expectations of the WIRE team. WIRE 
was to survey at least 170 square degrees to a lim- 
iting 12pm sensitivity of 1.9mJy, and smaller areas to 
deeper limits. Our source counts imply a surface den- 
sity of approximately 100 galaxies per square degree 
to the projected WIRE 12/im bright survey limit, us- 
ing the Pearson & Rowan- Robinson or Franceschini 
et al. 15/xm source count models to extrapolate. This 
is larger than projected source density of the WIRE 
team (q.v. figure 0, bottom-right panel, and @ left 
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panel, where the WIRE predictions were renormalised 
upwards by a factor of ~ 2 to match e.g. the observed 
IRAS and ELAIS counts). Although encouraging for 
the numbers of sources, it suggests that surveys much 
deeper than this will very rapidly be confusion lim- 
ited. This is in good agreement with the observations 
of Oliver et al. (1997). The composite counts from 
the various 15/xm surveys also imply that confusion 
will be irrelevant e.g. for NGST, but will be signif- 
icant for the MIPS imager on SIRTF even in short 
(~ 2000 second) exposures. The gains in detector size 
and sensitivity in future missions will offer very large 
improvements in wide-field survey efficiency: for ex- 
ample, a survey of a similar depth and areal coverage 
to ELAIS at 15/xm will be possible with MIPS on 
SIRTF in only ~ 10 hours. 

The upturn in the counts at faint fluxes continues 
in the Hubble Deep Field ISO counts, and is too large 
to be attributable to starburst K-corrections (e.g. 
Xu et al. 1998, Elbaz et al. 1998a,b). This excess is 
above any reasonable no-evolution predictions (figure 
^|), which all have shallower slopes than the evolving 
models independent of K-correction and world mod- 
els. This implies that the faintest sources in ELAIS 
are sampling a significantly cosmologically evolving 
mid-IR population. Given the strong evolution in 
both the comoving star formation rate claimed from 
optical-UV observations and in the quasar comoving 
number density, as well as the large fraction of AGN 
and starbursts expected in ELAIS from the counts 
models, we can reasonably expect that optical spec- 
troscopy of ELAIS and fainter samples will have a ma- 
jor impact on the study of dust-shrouded star forma- 
tion and quasar activity and their evolution. The pre- 
diction from the Pearson & Rowan- Robinson (1996) 
source count models is for of order a few objects as 
(intrinsically) luminous as IRAS F10214+4724 out to 
2 = 4 over the entire ELAIS areas. The MIR lumi- 
nosity function and luminosity density can be used 
to derive a comoving star formation rate without the 
difficult problems of extinction corrections that affect 
optical-UV estimates, and (for the ISOCAM LW3 fil- 
ter) is reasonably free of K-correction uncertainties 
(Xu et al. 1998, Serjeant 1999). The ELAIS survey is 
also highly sensitive to obscured quasars and will be 
an exceptional resource for active galaxy unification 
models. The ELAIS limits can also provide important 
constraints for sources detected at other wavelengths: 
for example, the 15/im, 90^im and 175/im limits at 
850/im source positions can provide robust redshift 
constraints (e.g. Hughes et al. 1998). 



5 CONCLUSIONS 

The extragalactic source counts agree extremely well 
with all evolving model predictions (Franceschini et 
al. 1994, Pearson & Rowan-Robinson 1996, Guider- 
doni et al. 1997, Xu et al. 1998) over seven orders 



of magnitude in 15pm flux density. The Pearson & 
Rowan- Robinson (1996) models can broadly repro- 
duce the 6.7/im extragalactic counts, but the obser- 
vations are in excess of the Franceschini et al. (1997) 
predicted counts at this wavelength using our prelim- 
inary 6.7/im flux calibration. All no-evolution mod- 
els are clearly excluded, and imply a cosmologically 
evolving population of obscured starbursts and/or ac- 
tive galaxies dominates below ~ 10 mjy at 15/xm, in- 
dependent of K-correction assumptions. Source con- 
fusion appears to have been underestimated in the 
WIRE and SIRTF missions, but will not impact sig- 
nificantly on the NGST. 
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